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1. INTRODUCTION 

Telecommunications allows transmissions of data over long distances with constraints of low power 
consumption and the low error rate with the help of error correction codes (ECC). Supported by the information 
theory and coding theory, ECC is playing important role and is highly required to detect and correct errors 
caused by channel distortion and/or noise corruption [1]. According to Shannon channel coding theorem [2], 
the probability of error can be made arbitrarily small or close to zero, if the transmission coding rate R is less 
than or equal to the channel capacity C’. Assuming that capacity C' is achievable, we can calculate the minimum 
signal-to-noise power ratio (SNR) to reach zero error or arbitrarily small error, called Shannon limit, if R < C 
is kept. If codes have close performance to the Shannon limit, the codes is said to have efficient transmissions, 
since for the given bit-error-rate (BER) level, the required SNR is smaller compared to the performance of 
uncoded case [2], [3]. 

Low density parity-check (LDPC) codes are a class of linear block codes in which the complexity 
of iterative decoding increases linearly with block-length [4]. As a result, LDPC codes have near-Shannon 
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limit performance on long block-length data transmission [5]. Gallager introduced regular LDPC codes with 
an equal distribution of ”1” in each row and column of the parity-check matrix H [6], whereas Luby developed 
irregular LDPC codes with an unequal distribution of ”1” in the matrix H. It has been proven that irregular 
LDPC codes perform better than regular LDPC codes [7], [8]. 

The idea of using LDPC-like codes for short block-length transmission, e.g., Internet-of-Things appli- 
cations has been proposed by [9] with Raptor codes based on LDPC and based on low density generator 
matrix (LDGM), of which the optimal the degree distribution is investigated in [11]. However, the applications 
of LDPC for large block-length has been investigated to in and for Digital Video Broadcasting Ter- 
restrial 2nd generation (DVBT2) and recently for the fifth telecommunication generation (5G) new radio and 
14], where Raptor-like structure is used to have rateless coding scheme capability. Recently, we revealed that 
irregular degree distribution is more beneficial since we have additional degree of freedom to construct better 
LDPC-like coding scheme as in [15]. In this paper, we extract those benefits of irregularity to design better 
codes for short block-length transmissions. 

The decoding process for LDPC codes uses a sum-product algorithm exchanging the log-likelihood 
ratio (LLR). This algorithm is represented by a Tanner graph involving variable nodes and check nodes. 
The check node receives the LLR from the variable node, which is subsequently returned to the variable node 
via a box-plus operation. This process is repeated until we obtain the desired outcome [18]. 

There are numerous problems in the data transmission process. Efficient data transmission is indicated 
by the low power consumption, but keep minimum errors [19]. Higher redundancy coding tends to perform 
better, but as the size of the data conveyed grows larger, so does the amount of power required [20]. Therefore, 
we introduce Stairs LDPC (SLDPC) codes, obtained from modification of the general simple LDPC for short 
block-length, to demonstrate the effectiveness of the equivalent structure derivation method. The coding struc- 
ture is expected to help source compression to increase transmission efficiency. We also do this to increase R, 
resulting in a more reliable transmission system. However, because this can affect the error correction perfor- 
mance of the proposed codes, it is important to investigate which structures have smaller gap to the Shannon 
limit. 

Therefore, this paper proposes equivalent structures for various compression structures to minimize 
encoding complexity. The equivalent structure is created by combining the two encoding procedures into a 
single, simpler step while keeping the transmitted bits unchanged. It is worth discussing here to investigate 
whether utilizing the equivalent structure improves SLDPC codes performance or not. The contributions of 
this paper are summarized as follows: 


1. The idea of stairs in the SLDPC codes is presented that the parity check is used to combine” the infor- 
mation bits such that higher coding rate is achieved. 


2. The derivation of complex graph is presented using cascaded generator matrix and parity check matrix 
such that the equivalent structure is obtained to minimize the complexity the encoder and receiver. We 
found that complex graph can be simplified into a simpler graph having the same performance. 


3. Validity of the derivation is confirmed using performance evaluation in terms of BER performance. We 
found that single parity check (SPC)-like codes are good for high coding rate. 


The rest of this paper is organized as follows. Section 2 presents the system model of the transmission, 
where additive Gaussian noise (AWGN) channel is the main channel used for the evaluation. This section also 
propose the SLDPC codes with 4 structures to be evaluated and compared with the 3 well-known existing codes 
structures. Section 3 evaluates the proposed codes performances followed by conclusions in Section 4. 


2. METHOD 

In this section, we separate the discussion of methods into two subsections. Subsection[2.1]] presents 
system model of the transmitter and receiver structure. Subsection [2.2] presents method to derive equivalent 
coding structures of elements to construct LDPC codes. 
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Figure 1. Transmitter and receiver structures of SLDPC codes. 


2.1. System Model 

A set of tests are presented in this paper to examine the error-correcting capability of SLDPC codes 
throughout data transmission. To maximize the coding rate, we introduce a concept of a compression structure 
followed by the equivalent structure. Figure [iJillustrates the system model considered in this paper. 

The information bits u are random binary bits encoded in channel coding block Cg with a block length 
and various coding rates. The encoded bits x are then modulated by / to obtain s prior to the transmission by 
the antenna. At the receiver, the signal y experiences the channel h and corrupted by noise n resulting in 


een (1) 


The LLR of the channel 9 
Leh = sr tty} (2) 


is then demapped to obtain Le,jy with R{-} being the operation to take the real part of y. The Ley is then 
used by the SLDPC decoder Dg resulting L;. Finally, hard decision is performed to obtain w. 

We use binary phase shift keying (BPSK) modulation since it has a simple modulation mechanism 
and a longer transmission range. This process changes the bit ’0” to 1” and bit ”+1” to ”—1”. The result of 
BPSK modulation is given by 


-lz= 


The modulation results s are then transmitted across the transmission channel. 

We use additive white Gaussian noise (AWGN) channel as the transmission channels, however, exten- 
sion to fading channel is rather straightforward by modifying (i). The terms of additive refers to noise that is 
simply superimposed or added to a signal without any kind of multiplication mechanism and is the same 
regardless the frequency band inside the used bandwidth. Theoretically, the bit-error-rate (BER) of uncoded 
BPSK under the AWGN channels can be expressed as 


1 | E 1 / 1 
P..AWGN => gente ( 2 => gente ( a 5 (4) 


with fn being the energy bit per noise, which is used to confirm the validity of the computer simulation and 
the baccline of the performance comparison. The variance is given by [22] 


o = 1Q~7[4B]/10 (5) 


with 7 being the SNR in dB. Since we use BPSK modulation, Le,jz equals to L-y. Interested readers on 
mapping other than BPSK can refer, for example, the demapper for bit-interleaved coded modulation (BICM). 
The L, 7 1s then used as the input of soft decoding process to return back information t from the received 
signal. 

The sum-product algorithm is one of the optimal iterative decoding techniques in the logarithmic 
domain based on soft-decision. The ’sum” represents operations on variable node (VND), which is the sum- 
mation over all incoming messages. On the other hand, ’ product” represents operations on check node (CND), 
which is expressed using a ’box-plus” operation, that perform a special multiplication over all the incoming 
messages. Figure [2] shows the operations for VND and CND, of which the outgoing message is Le vnp and 
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LeVND 


(a) 
Figure 2. Principle of operations for (a) VND and (b) CND for equations in (oh and (8). 


Le,onp> respectively. This figure also shows that L,;, is only connected to the VND, while the CND only 
depends on the Ly vnp. 

The operation on VND is divided into two parts, i.e., (i) extrinsic LLR and (ii) a posteriori LLR. The 
extrinsic LLR for the 7-th edge outgoing from the VND is expressed as 


dy 
Le,vnp,i = Len + S- LaNND,j- (6) 
j=l j#i 
The symbol Le,vnp,; is the j-th a priori LLR coming to the VND from CND with exception of 7 # i. The a 


posteriori LLR is calculated in the end of decoding after the number of iterations is assumed to be enough and 
is expressed as 


d 


Lynd = Len + S- LaND,j: (7) 
j=l 


where we omit the index 7, since for the final process, we need only a single LLR for each desired symbol. 
The box-plus operation at CND is expressed as 
de 


Lecnpi= >, WLa,onv,j; (8) 
j=l ii 


where the box-pus arithmetic operation is given by 


L, & Lz = sign (L, - Lz) min (|L£4], |L2]) (9) 


for the case of two edges. Referring to (9), the box-plus operation is performed by multiplying the sign of Ly 

and Lo, then multiplying it by the smallest absolute value of L, and L» [23]. This equation indicates that zero 

LLR coming to the CND should be avoided since all results are becoming zero causing no benefit in decoding. 

This is one of the reason that LT codes require more protections, for example, by using accumulator as in (15). 
The outcomes of L,, vnp are then returned to bits using hard decision of 

A 0, Lip .WND >0 

ai { 1, Lp vnp < 0, a) 


where %; is i-th bit of final decoding. In this paper, we use hard decision based on LLR as follows. If Lp vnp 
is greater than 0, w is 0. On the other hand, if Lp vnp is less than 0, % is 1. This is to be consistent with the 
mapping rule in (3). 

However, because the focus of this paper is on the design of SLDPC codes with high coding rate, we 
start the design from the encoder graphically via a generator matrix G before converting it to matrix H using 


G=[|-I, | P], H=[P? | Tn-x] (11) 
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(a) 
Figure 3. Illustrations of (a) the stairs principle and (b) their realizations for data compression by CND. 
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Figure 4. The graphical structure of SLDPC1 codes. 


where G is k x n generator matrix, I, is k x k identity matrix, and P is k x (n — k) designed matrix. This 
technique is similar to the design of LDGM codes [9], [10], however, in this paper, we seek the high coding rate 
and small block-length. The scalability of the codes can then be extended by using protograph-based extension 
technique [24]. 


2.2. Deriving the Equivalent Coding Structure 

In this section, we present a method to derive equivalent coding structures of elements to construct 
LDPC codes using SLDPC as a toy example. The principle idea of SLDPC codes is presented and shown 
in Figure [3] Figure Be) illustrates the stairs principle representing a compression encoding to obtain higher 
coding rates to help transmissions of source compressed data. Figure [3{b) realizes how the CND can be used 
to compress the data such that higher coding rate is achievable. 

We present SLDPC codes having 4 structures compared to the existing well-known coding schemes 
as shown in Figures [4}{10} which are then summarized in Table[1] All structures are having unique properties 
and is intended for data compression such that rate is kept high. Please note that in this paper, we start the 
discussion on the structure from the generator matrix G with the graphical viewpoint prior to the conversion to 
parity check matrix H. 


2.2.1. SLDPC1 Codes Structure 
Structure 1 of SLDPC codes, called SLDPC1 codes, shown in Figure |4} is the initial structure of the 
proposed SLDPC codes. Based on the structure of encoding processes, we split the generator matrix into 


1041 1 1 0 1 
ieee 1 ile (12) 


where Gy, is for the left part graph and G1, is for the right graph at the transmitter side. We are now ready to 
convert matrices G,; and Gj,. into matrices Hj; and H},, respectively, based on (ip. The matrices H,; and 
H,,. are therefore becoming 


1 11 0 


Hu =|; 001 


|. Hi,=[(1 1 1]. (13) 


Since performing two decoding processes in are considered to have higher decoding complexity, 
we propose to simplify it using the concept of Cascaded LDPC codes [25], such that we only have both a single 
equivalent Gyeq and Hyeg. 
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We start from the equivalent matrix Hj, since it is easy to construct as 
1 110 0 
Hi = is ale 1001 of, (14) 
00111 


where the O, has the size of k x 1, while 0, has 1 x & adjusting the size of Hy; and Hj, where k = 2. 

The matrix Hj. represents VND and CND of the SLDPC1 codes. The number of VND and CND is 
equal to the number of columns & and rows n, respectively. Refering to (14), SLDPC1 codes are in the class of 
irregular LDPC codes due to unequal distribution of ”1” in the matrix Hy eq. 

Similar to the case of parity check matrix case, the encoding generator matrix is also simplified into 
the equivalent matrix G_,, which is obtained from Gy; and G,,. The matrix G_g can be obtained based on the 
cascaded principle as 


(15) 


G 0 
Gi cascaded = | _ a _ 


0, Gi, 


ooo F 
ocooro 
OrRrRH 


Pil p2 


where the O, has the size of k x 1, while 0, has 2 x k following the number of row of G1; and number of 
column of G,,, where k = 2. 
The input length is only 2 bits, for example b; and bg. The parity bits generated by Gy, are, therefore, 


pi = b @be, (16) 
p2 = by. (17) 


We then further adjust the matrix by locating the last parity check bit of G, in the last column of Gjeq such 
that the input length is kept 2 bits. The matrix G1, indicates that the last generated parity check is 


p3 = pi © p2 = (b1 G ba) @ by = by. (18) 


Therefore, the last parity check in the matrix G_, is [0; 1)” selecting only the bit bj. The final equivalent 
generator matrix is 


1 01 1 0 
Cred = 0 110 | oe 
Based on (19), the parity check matrix Hj. for SLDPC1 codes is 
1 1100 
Mie = ]1 0 0 1 OF}. (20) 
0 1 0 01 


The last step is to confirm, whether the matrices in (20) and (14) is the same or not, since they are 
obtained from two different ways. Hye, in is obtained from cascaded Hy, and H1,, while Hi¢q in is 
from the cascaded G,, and Gj,. Here, with Gauss-Jordan elimination by XOR-ing: (i) the first and the third 
rows followed by (ii) the second and the third rows, we obtain 


1 11 00 1 1 1 0 0 1 11 0 0 
Meg = {1 0 0 1 OJ =]1 0 0 1 OJ =]1 0 O 1 OF} = Meg (21) 
0 10 0 1 101 0 1 00111 


confirming the equality of Hyeq° = H1eq and validity of the equivalent matrix. The coding rate of SLPDC1 


codes is Ry = 2, since 2 bits are encoded into 5 bits as shown in Figure|4] 
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Na A A 
h Uy ug 


uy 


ug 


Figure 6. The graphical structure of SLDPC3 codes. 


2.2.2. SLDPC2 Codes Structure 


To obtain many possible high coding rate, we also propose the second structure of SLDPC codes, 
called SLDPC2 codes, as shown in Figure 5] with the coding rate of Rg = : = 5. This SLDPC2 codes are 
obtained from SLDPC1 codes by removing the last parity check of G1,. The generator matrix is 


1. 0.4. 4 
G =|) a 6] = Gun (22) 


Please note that the VND in Figure |5]can be ignored since the degree is two, which is similar to a line. The 
parity check matrix for SLDPC2 codes is 


= Mir. (23) 


It is interesting to observe here, whether the SLPDPC1 codes with additional parity check ps3 has better perfor- 
mance or not compared to this SLDPC2 codes. 


2.2.3. SLDPC3 Codes Structure 


We also propose the third SLDPC codes, called SLPDC3 codes, by modifying the SLDPC1 codes 


with removal the link of VND but keep the parity check CND as shown in Figure [6] The channel coding rate 
of SLDPC3 codes is Rg = z = $3 of which the generator matrix is the same as Gieg in except in the 


forth column, which should be removed, indicated by fer in the 4-th column, due to the deletion of the link 
connecting VND and antenna. Therefore, the equivalent matrix is 


101 +4 0 1 01 0 
Grea = [5 1106 |= (5 11 i ae 
The parity check matrix is, therefore, 
1 11 £0 
Hoea = |4 i a ik (25) 


Please note that the matrix G3-, and Gg is similar except at the last column suggesting their performances 
comparison is of our interest. 
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Figure 7. The graphical structure of SLDPC4 codes. 
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Figure 8. The graphical structure of Repetition codes. 


2.2.4. SLDPC4 Codes Structure 

To further reach the high coding rate, we further modify the SLPDC3 codes by removing the link 
Connection VND to the antenna, called SLDPC4 codes, such that the channel coding rate is Ry = 3 as shown 
in Figure This SLDPC4 codes structure has the highest channel coding rate R among all the proposed 
SLDPC codes. The generator matrix is obtained by further deleting the third column of G3¢,, indicated by 


+4" in the 3-rd column, as 


1 0 + 0 1 0 0 
Gaeq = F 1 4 | _ i 1 | oo 
The equivalent parity check matrix is, therefore, given by 
Hig = (0 1 1s (27) 


It is important to note here that the evaluation on the performance of SLDPC4 codes is important since H4eq 
has zero element resulted, which is in general useless. 


2.2.5. Repetition Codes Structure 

To make our investigation complete, we also provide comparison with the existing famous Repetition 
codes, since the Repetition codes are simple and have easy adjustable coding rate. The structure of Repetition 
codes with rate R5 = 3 is shown in Figure |8] with the generator and parity check matrices are, respectively, 
given by 


(28) 


Gy= [1 1 a] and rig=[t 2 9). 


2.2.6. LDGM Codes Structure 

In this paper, we also investigate a comparison of the proposed SLDPC codes with the existing well- 
known low density generator matrix (LDGM) codes. To keep comparable simple structure, we use LDGM with 
Re = 2 as shown in Figure|9| The generator and parity check matrices of the considered LDGM codes are, 


5 
respectively, expressed as 


1110 0 
Ge=|j 1 1 1 | md He=|0 10 1 0 (29) 
10001 
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Table 1. Summary of the proposed SLDPC codes and the existing well-known codes. 


No | Codes Matrix G Matrix H Rate R 
1110 0 
Hig =|1 0.0 1 6 
1 1 i 11 
1 | The proposed SLDPC1 | Gicg = : ’ aa 2 
01101 tk 2 Bo 
Hiew=|1 00 1 0 
01001 
130 1 1 1 110 
2 | Th d SLDPC2 | G, = H, = 2 
lee : F 4 \ : ; 0 0 4 
1 1 0 1 110 
3 | Th d SLDPC3 | G3., = H3, 7 
€ propose 3eq 0 3eq f l 0 i Zq 
1 0 0 2 
4 | The proposed SLDPC4 | Gieg= | 1 Hieq= [0 1 1] a 
1 41-6 
5 | Repetiti Gs=[1 1 1] H; = i 
epetition 5 5 ; 0 i 3 
it de Ura 1 i206 
6 | LDGM cs| | H;=|0 10 1 0 Z 
01110 5 
10001 
ie | 
7 | SPC G;= H =(1 1 1] 2 
: f 1 i r 3 


2.2.7. SPC Codes Structure 


Since the well-known single parity check (SPC) codes are having high coding rate of R = ei with 
k; being the information bits, we also investigate the performances of the proposed SLDPC codes compared 
to the SPC codes. Figure[10|shows the SPC codes structure used in this paper, where the channel coding rate 
is R7 = 2, since we select k = 2 bits. The generator and parity check matrices of the (2,3) SPC codes are, 


respectively, given by 


t- Oa 
G=|j i i and H;=[1 1 1). (30) 


We summarize all coding schemes and their corresponding channel coding in Table [I] with their cor- 
responding generator and parity check matrices G and H, such that their differences are also comparable and 
visible. 
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Figure 9. The graphical structure of LDGM codes. Figure 10. The graphical structure of SPC codes. 
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Figure 11. BER performances of SLDPC codes under AWGN channels. 


3. | RESULTS AND DISCUSSION 

In this section, we present performance evaluation in terms of BER under AWGN channels. We have 
conducted a series of computer simulations to figure out the performances of the proposed SLDPC codes and 
find the best SLPDC coding structure among 4 possible structures. We use signal-to-noise power ratio (SNR) 
ranging from —7 to 10 dB under the AWGN channels and the corresponding coding rates R in Table[I] Asa 
baseline comparison, we use the theoretical BER performances of AWGN channel to confirm the validity of 
our simulations and the performances comparison to the proposed SLPDC codes. We present the simulations 
result into two figures of BER performances to make the BER curves visible and comparable. We also plot the 


Shannon limit of 
‘tim = 10 - logy)(2” — 1) (31) 


with yLim being the SNR minimum to be theoretically achieved with the given channel coding rate R from 


channel capacity 
C= B-log,(1+7) (32) 


with B being the bandwidth to fairly evaluate the efficiency of each coding structure. The limit is said to be 
achieved when C' = R for arbitrarily small BER. Figure [II] depicts the BER performances of SLDPC1 with 
Hy.g, SLDPC1 with H,.4, Repetition, and LDGM codes, while BER performances of SLDPC2, SLDPC3, 
SLDPC4, and SPC codes are shown in Figure [12]completed by the Shannon limits for the corresponding rates. 

We performed soft iterative decoding for all coding structures. The iteration process is carried out 
by repeating the sum-product algorithm until a steady LLR value is reached. Since the other structures had 
reached a stable LLR value by the time of the first iteration, we only repeated the iteration process for SLDPC1 
and LDGM codes. Figure [IT|shows BER performances evaluated in terms of SNR(dB), where the Repetition 
codes has gap of 9.35 dB to the Shannon limit, while SLDPC1 codes with Hj¢q’, LDGM codes, and SLDPC1 
codes with H1¢q have the same gap of 8.45 dB. All gaps are measured at BER level of 10~*. The results show 
that different bipartite graphs of SLDPC1 codes and LDGM codes can provide the same BER performances 
indicating that the proposed method of equivalent structure derivation is correct. 

Figure[i2]shows that SPC codes have the best performance indicated by the smallest gap of 8.01 dB. 
SLDPC2 and SLDPC3 codes have better BER performance compared to SPC but have less efficiency. Figure[12| 
also shows that SLPDC2 and SLDPC3 codes have the same performances because the generator and parity 
check matrices is the same. The SLDPC4 codes have high rate but the performance is the worst with gap of 
10.31 dB. This is because H4eq in has zero element for a single parity check meaning that LLR to that 
VND is useless. 
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Figure 12. BER performances of SLDPC codes under AWGN channels. 


4. CONCLUSION 

In this paper, we have proposed a method to derive equivalent coding structures of elements using 
SLPDC codes that can be extended to help the construction of larger LDPC codes for short block-length trans- 
mission applications. We used the cascaded principle to construct the equivalent generator and parity-check 
matrix that is also beneficial for minimizing the computational complexity of the encoder and decoder. The 
equivalent structure of G can be obtained from the cascading of two (or more) matrices by modifying the sec- 
ond matrix to have identity matrix with size of parity check length and shifting to the right k bits. On the other 
hand, the equivalent structure of H can be obtained by just cascading directly the second matrix and shifting to 
right k bits. We have evaluated the validity of the proposed method by confirming the equality in performances 
of the SLDPC codes and other well-known coding scheme under the AWGN channel. According to the sim- 
ulation results, the proposed method is valid since SLDPC1 codes have the same BER performance with that 
of LDGM codes. Although some bipartite graphs are different to each other, they can share the same equiva- 
lent matrix and provide the same performances indicating that the equivalent structure derivation is important. 
Among all the channel coding structures evaluated in this paper, we found that SPC codes are the most efficient 
codes (for the given block-length) since SPC codes have the smallest gap to the Shannon limit. The results in 
this paper are expected to provide new insight for the designing simple channel coding of short block-length 
LDPC codes having high coding rate for future less power consumption applications. 
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